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Abstract 

(N 

The Card-Cyclic-to-Random shuffle on n cards is defined as follows: at time t remove the card 
with label t mod n and randomly reinsert it back into the deck. Pinsky [3] introduced this 
shuffle and asked how many steps are needed to mix the deck. He showed n steps do not suffice, 
i— h Here we show that the mixing time is on the order of Q(nlogn). 

P 
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, — , 1 Introduction 

PLh In many Markov chains, such as Glauber Dynamics for the Ising model, the state space is a set of 

configurations, and at each step a location is chosen and updated. An important general question 
about such chains is what happens when we move from the world of random updates, where at each 
step a location is chosen at random and updated, to systematic scan, when the updates are done in 
a more deterministic fashion (see e.g., [1]). On the one hand, systematic scan is "less random", so 
one might expect that the mixing time is larger. On the other hand, systematic scan can update 
J> n sites in n steps, whereas with random updates n log n steps are required by the coupon collector 

problem, so one might expect systematic scan to have a smaller mixing time. 

This question has been investigated in the context of the random transpositions shuffle. In this 
shuffle, at each step a pair of cards is chosen uniformly at random and interchanged. In a classical 
result of Diaconis and Shahshahani [2j, the mixing time of the random transposition shuffle is 
shown to be asymptotically ^nlogn. Mironov [7], Mossel, Peres and Sinclair [8j, and Saloff-Coste 
and Zuniga [llj analyzed the Cyclic-to- Random shuffle, which is a systematic scan version of the 
random transposition shuffle: at step t the card in position t mod n is interchanged with a randomly 
chosen card. They found that the mixing time for this chain is still on the order of nlogn. 

In the present paper, we study a systematic scan version of the Random-to-Random insertion 
shuffle. In the Random-to-Random insertion shuffle, at each step a card is chosen uniformly at 
random and then inserted in a uniform random position. It was shown in [3], |12j and [10] that 
the mixing time of this shuffle is on the order of n log n. Pinsky [9] introduced the following model, 
called the Card-Cyclic-to-Random shuffle: at time t remove the card with the label t mod n and 
insert it in a uniform random position. It is not obvious that the mixing time is greater than n: 
after n steps the location of each card has been randomized, so one might expect the whole deck 
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to be close to uniform at time n. However, Pinsky showed that the mixing time is indeed greater 
than n, since the total variation distance to stationarity at this time converges to 1 as n goes 
to infinity. We show that in fact the mixing time is on the order of nlogn. To prove the lower 
bound we introduce the concept of a barrier between two parts of the deck that moves along with 
the cards as the shuffling is performed. Then we show that the trajectory of this barrier can be 
well-approximated by a deterministic function / satisfying 

/(*) = f f{s)ds, (1) 

Jx-l 

and we relate the mixing rate of the chain to the rate at which / converges to a constant. To prove 
the upper bound, we use the path coupling method of Bubley and Dyer pp. 



2 Statement of main results 

Let Xt be a Markov chain on a finite state space V that converges to the uniform distribution. For 
probability measures \i and v on V, define the total variation distance — z/|| = \ Ylx&v ~~ 
and define the e-mixing time 

imix(e) = min{t : \\F x (X t = •) - U\\ < e for all x € V} , 

where IA denotes the uniform distribution on V . 

Recall that in the Card-Cyclic-to-Random shuffle, at time t we remove the card with label 
t mod n and then reinsert it into a uniform random location. 

Define a round to be n consecutive such shuffles. Note that the Markov chain that performs 
a round of the Card-Cyclic-to-Random shuffle at each step is time-homogeneous with a doubly- 
stochastic transition matrix, irreducible and aperiodic, hence converges to the uniform stationary 
distribution. It follows that the Card-Cyclic-to-Random shuffle converges to uniform as well. Our 
main results show that the mixing time is on the order of logn rounds. 

Theorem 1 There exists cq such that for any c < cq and < e < 1, when n is sufficiently large, 
we have 

imix(e) > cnlogn. 

Here c$ = 2 +2a w ^ iere a ^ s ^ e smallest positive solution of equations b = e a sin b and a = e a cos 6—1. 
Numerically c$ = 0.161875162.... 

Theorem 2 For any e > and n > 4, we have 

imix(e) < C(n log n - 2n log e), 

where 

C = —. r = 6.58664655... 

log 2 — log(e — 1) 

Remark. Theorem [T] and Theorem [2] together establish that the Card-Cyclic-to-Random shuffle 
has a pre-cutoff in total variation distance. It is an interesting open problem to determine if cutoff 
occurs in this shuffle. For reference on cutoff phenomenon, see [BJ Chapter 18]. 
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Remark. One can also consider a simpler shuffle where at time t, the card in position t mod 
n is removed and inserted in a uniform random position. Call this Position-Cyclic-to-Random 
insertion. The time-reversal of this chain can be obtained by at time t, picking a uniform card 
and inserting it to location t mod n. Considering the length of the longest increasing subsequence 
shows O(nlogn) steps are needed to mix(Ross Pinsky, personal communication). A matching 
upper-bound of 0(n log n) follows from the work of Saloff-Coste and Zuniga. See Theorem 
4.8]. 

Theorems [I] and [2] will be proved in Sections [3] and |4j respectively. 

3 Lower bound 
3.1 The barrier 

The key idea for the lower bound is to imagine a barrier between two parts of the deck, that moves 
along with the cards as the shuffling is performed. If a card is inserted into the gap that the barrier 
occupies, we use the convention that the card is inserted on the same side of the barrier as it was 
in the previous step. We illustrate this with the following example. Suppose there is a deck of 8 
cards with a barrier between cards 3 and 5. In the next step, card 7 is inserted between cards 3 
and 5. 

2 1 3 | 5 4 6 8 7 
2 1 3(7 5 4 6 8 

Let {oi}^ be a Card-Cyclic-to- Random shuffle. We think of o~t(i) as the position of card i at time 
t, where the positions range from 1 at the left to n at the right. Define the position of the barrier 
as the position of the card immediately to its left, and throughout the present chapter, let Bt be 
the position of the barrier at time t. Use the convention that Bt = if at time t the barrier is to 
the left of all cards. We will call the pair process (at, B t ) the auxiliary process. 

Note that the conditional probability that the card at time t is inserted to the left of the barrier, 
given Bt, is r-B*- Since at any time t > n, every card has been moved exactly once in the previous 
n steps, we have 

n 

Bt = 1 (th e card moved at time t — i is inserted to the left of barrier) , 
i=i 

and hence 

n 

E(B t ) = - VlfBi-j). (2) 
n 

i=i 

Define g(t) = K(^Bt). Then g satisfies the following moving average condition: 

1 n 

= -!>(*-*')> ( 3 ) 

i=i 

for t > n. We shall approximate g(t) by f(t/n), where / : R — > [0, 1] is a continuous function 
satisfying (TTj) . Our first lemma gives an example of such a function. 
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Figure 1: Graph of f(s) 



Lemma 3 There exists a > and b > 2ir such that f(x) = \ + \e ax sin(6x) satisfies 

f'(x) = f(x)-f(x-l). 



(4) 



Moreover, 



f(x) = / f(s)ds, 

lx-l 



(5) 



for all x. 



Proof: Since properties Q and §5§ are preserved under shifting and scaling, it is enough to show 
that they apply to h(x) = e~ ax sin(6x), for suitable a and b. 

First, we show that for suitable choice of a and b we have h'(x) = h(x) — h{x — 1). By the 
product rule, 

h'(x) = -ae~ ax sm(bx) + be~ ax cos(te), (6) 

and a calculation shows that 

h(x) - h(x - 1) = (1 - e a cos6)e- ax sin(6x) + (e a smb)e~ ax cos(6x). (7) 

The quantities ^ and ([7]) are equal if b = e a sin b and —a = 1 — e a cos b. Solving for a in the first 
equation gives 

, b 
a = log ■ 



sin b' 



and substituting this into the second one gives 



sin b b cos b 

log—— = 1 — -. 

o sin b 
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By the intermediate value theorem, this equation has a solution with b in the interval [27r+f, 2tt+^], 
since when b = 2tt + ? the right-hand side is smaller than the left-hand side, but when b = 2ir + | 
the right-hand side is larger. Furthermore, since sin b < b when b > 0, we have a = log -Ag > 0. 
(Numerical approximation gives the solution as b = 7.4615... and a = 2.0888....) 

Next we claim that since h'(x) = h(x) — h(x — 1), we must have h{x) = J h(s)ds. To see 
this, define q(x) = f _, h(s)ds and note that = h'(x). This implies that — g(x) = C for 
a constant C. But since o > 0, we have h(x) — >• as x — > oo. Consequently — )• as x — > oo by 
the definition of q, and so C = 0. □ 

Recall that = K(-^Bt), where Bt is the position of the barrier at time t. A key part of our 
proof will be to show that g closely follows the continuous function / of Lemma [3j However, in 
order for this to be the case we must start with a permutation chosen from a certain probability 
distribution. It is most convenient to describe this starting permutation as being generated in the 
first n time steps, which we call the startup round. In the startup round, we begin with only a 
barrier. At time t, for 1 < t < n, we put card t to the left of the barrier with probability /(^)- 
The location among the already existing cards in the left (right) side of the barrier is arbitrary. We 
must modify the definition of g to handle the startup round. Define g : {1, 2, ...}—)• R by 



g(t) 



fO if 1 < t < n; 
E(-^Bt) otherwise. 



Thus g satisfies the moving average condition, and, because of the insertion probabilities used in 
the startup round, g matches / for the first n steps. (That is, g(-) = /(-) on {1, • • • ,n}.) As we 
show below, this is enough to ensure that g is well-approximated by / for a number of rounds on 
the order of log n. 

Lemma 4 There exists a constant C > such that 

\9(t)-m\<f n e 2{t+1)/n , 

for all t > 0. 

Proof: First, note that if t > n then 



g(t+l)-g(t) = IJ^t + l-O-IJ^t-i) 

i=l i=l 

= l(g(t)-g{t-n)). 

Rearranging terms gives 

g(t + l) = (l + ±)g(t)-±g(t-n). (8) 

Recall that f(x) = \ + \e~ ax sin(6x) and a > 0. Some calculus shows that the second derivative of 
/ is uniformly bounded on [0,oo). Hence 
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where the first line follows from Taylor's theorem and the second line follows from Lemma |3j 
Rearranging terms gives 



/(^) = (i + i)/(^-i/(^) + o(^)- 

Combining Q and ^ and using the triangle inequality gives 

\g(t + 1) - fm\ < (i + £) - + J |*(t - n) - /(*=*)! + §, 

for a universal constant C. We claim that for all t we have 

M*)-/(£)|<£E(i + S)*. 

i=0 



(9) 
(10) 

(11) 



We prove this by induction. For the base case, note that g(t) = f(^) for t = 1, . . . ,rt. Now if we 
suppose that ([Til holds for 1, . . . , t, 

be bounded by^EUoCl + s)*- H 



then the two absolute values on the right-hand side of ( 10 ) can 
Hence 



\g(t + i)-/( 



n ' 



< (1 + 
t+1 



1 + 



2 \i 



+ C 



i=0 



2 \i 



i=0 



which verifies ( 11 ) for t+1. To finish the proof of the lemma, note that 



c 

n 2 



(1 + -) 



2\t+l 



i=0 



^ 2^ e 



CJ(+l)/n 



3.2 Deviation estimates 

In the previous subsection we proved that the expected barrier location is well-approximated by a 
continuous function. In the present subsection we show that the barrier stays reasonably close to 
its expectation with high probability when the number of rounds is on the order of log n. 

Define a configuration as a pair (a, b), where a is a permutation and b is a barrier location. (Thus 
the state space of the auxiliary process is the set of all configurations.) We define the insertion 
distance between two configurations as the minimum number of cards we would need to remove and 
re-insert to get from one configuration to the other. For example the insertion distance between 
the two configurations below is 2. (Move cards 4 and 7.) 

2 1 4 3 | 5 6 8 7 
2 1 3 7 | 5 4 6 8 

Lemma 5 Let (al,Bj) and (cr|,i?|) be auxiliary processes, and define a\ = (o~\ : B\) for i = 1,2. 
Let d be the insertion distance between ctq and ctq. Then 

pEBi-Eflfl <d(l + I)*. 
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Proof: There is a natural coupling of o~\ and a 2 that we call label coupling. In label coupling, at 
time t we choose a label X uniformly at random. If X = t mod n, then we move card t mod n to 
the leftmost position in both processes. Otherwise, we insert card t mod n to the right of the card 
with label X in both processes. 

Suppose that A = {a\, . . . , a^} is a minimal set of cards that can be moved to get from <Tq 
to <Tq. Note that under the label coupling, only in the case when we move a card not in A can 
the insertion distance be increased. In such moves, if the card is put to the right of a card in A, 
the insertion distance increases by 1 and otherwise it stays the same. Thus the expected insertion 
distance after one step is at most 

/, -,sd n - d / 1 

(d+l)- + d =d 1 + - 

n n \ n 

Iterating this argument shows that the expected insertion distance after t steps is at most d (l + ^) . 
The lemma follows from this, since the barrier can move by at most one position with each re- 
insertion. □ 

We are now ready to state the main lemma of this subsection. 

Lemma 6 Let (at, Bt) be an auxiliary process. Fix c > and suppose T satisfies n < T < cnlogn. 
Then for any x > we have 

P(| S-Bt - 9(T)\ >x)< 2exp (-xV~ 2c ) . 

Proof: Fix T with n < T < cnlogn. Since g(T) = -K(Bt), it is enough to show that for any 
x > we have 

F(\B T - E(B T )\ >x)< 2exp {-x 2 ^^ 2 ^ 
Let Tt be the sigma- field generated by the process up to time t, and consider the Doob martingale 

M t :=E(B T \T t ). 

Applying Lemma [5] to the case of two configurations that differ by one insertion gives 

\M t -M t -i\ < (1 + ^) T ~*, 
for t with 1 < t < T. Thus the Azuma-Hoeffding bound gives 

F(\B T -E(B T )\ > x) = F(\M T -E(M T )\>x) 

K - 2 -te)' (12) 

where bt = (l + h) ■ Let r = (l + -) . The sum in (12) can be written as 

T-l T . 

r _i 

8=0 

< §r T , (13) 



since r — l = ^ + ^2>^. Since T < cnlogn, the quantity (13) is at most 



Substituting this into (12) yields the lemma. 
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3.3 Proof of the Lower bound 



Recall that /(g) = \ + \e~ as sin(6s), for some a = 2.0888... and b = 7.4615.... The rough idea for 
the lower bound is as follows. Note that if c is sufficiently small and s < clogn, then the fluctuation 
of f(s) between s and s + 1 is of higher order than n^ 1 / 2 . Thus in the corresponding round of the 
Card-Cyclic-to- Random shuffle, there will be an interval of cards where the probability of inserting 
to the left of the barrier is detectably high. Before we give the proof, we recall Hoeffding's bounds 
in 0. 

Theorem 7 Let Xi, . . . , X/- be samples from a population ofO's and 1 's, and let p = K(X\) be the 
proportion of 1 's in the population. Then 

F \J2 Xi-kp>o)j < e~ 2a2/k . (14) 



The bound ( 14 ) applies whether the sampling is done with or without replacement. 
Proof of Theorem [1} Let c > be small enough so that 

C< 2TYa < 15 > 

Fix T with n < T < cn log n and let x = T/n. Suppose that sin(6x) < 0. The case sin(frx) > is 
similar. Since b > 2ir, there exist x\,X2 with x — 1 < x± < X2 < x, such that 

bx\ = 2-Kk + 7r/4, and 

bx 2 = 2nk + 37r/4, 
for an integer k. Note that for s £ [x\, xz] we have 

f(s) > \+Pe~ as 

> \+pn~ ac , (16) 

where f3 = \ sin(7r/4). The second inequality holds because x < clogn. 

Let A be the event that \±B t - f(t/n)\ < fn" ac for t with T - n < t < T. Note that 
since T < cn log n, substituting T into the upper bound of Lemma [4] implies that if t < T then 
\g(t) — f(t/n)\ < Bn 2 ^ 1 , for a constant B > 0. Since 2c — 1 < —ac by ( |l5[ ), for sufficiently large n 
we have 

Bn 2c - 1 < § n~ ac , 
and hence \g(t) — f{t/n)\ < f n~ ac for t < T. Hence 

F(A C ) < P(\±B t - g(t)\ > f n~ ac for some t with T - n < t < T) 
< 2n exp ( 



2nexp(-gn 1 - 2c ^ 1 )), (17) 



64' 

where the second inequality follows from Lemma [6] and a union bound. Since 1 — 2c(a + 1) > by 
(15), the quantity (17), and hence F(A C ), converges to as n — > oo. 
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Let I = {t mod n : nx\ < t < 11x2} and m = Since X2 — x\ = vr/26, there is a constant 
A such that m > An for sufficiently large n. Let N be the number of cards in / (that is, cards 
whose label is in /) placed to the left of the barrier between times nx\ and nxi- Then ./V is also the 



number of cards from / to the left of the barrier at time T. By (16), on the event A the insertion 
probabilities ^ are bounded below by \ + ^n _ac for t with nx\ < t < nx2- Hence the conditional 
distribution of N given A stochastically dominates the Binomial(m, \ + 
Hoeffding's bounds give 



n 



distribution. Thus 



N < f + fmn" ac | A) < exp 



1-2 ( |fJMJ~" r 



in 



< exp 



-An 



l-2oc 



(18) 



where the second line follows from the fact that m > An. Since 1 — 2ac > by (15), the quantity 



( 18) converges to as n — > 00. 



Now let Y be the number of cards in / having position less than § + ^n 1 ac at time T. Since 
/ (?) < |) w e have -Br < § + \n L ~ ac on the event ^4, and hence 



P(Y < f + f mn~ ac ) < P(A^ < f + fmn- ac |^) +P(A C ), 



(19) 



which converges to as n — >■ 00. 

To complete the proof, let 1^ be the number of cards in / whose position is less than ^ + ^n 1_ac 
in a uniform random permutation. 

Hoeffding's bounds imply that 



f + f mn" ac ) 



< exp 



2 ( fmn" ac 



< exp 



(20) 



for sufficiently large n. Since 1 



this with (19), we conclude that t 



2ac > 0, the quantity (20) converges to as n 
imxi e ) — cn log n for large enough n. 



00. 



Combining 

□ 



4 Upper bound 

We use the path coupling technique introduced by Bubley and Dyer [I]. Let S n be the permutation 
group and G = (S n ,E), where an edge exists between two permutations if and only if they differ 
by an adjacent transposition. The path metric on G is defined by 

p(x, y) = min{ length of r] : 7] is a path from x to y}. 

Define 

diam(G) = sup p(x,y). 

x,y 

The following theorem is from [IJ. See also [6J Chapter 14]. 
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Theorem 8 Suppose that there exists a > such that for each edge {x, y} in G there exists a 
coupling (Xi,Y±) of the distributions P(x, •) andP(y,-) such that 

E x>yP (X 1 ,Y 1 )<p(x,y)e- a . 

Then 

-logg + log(diam(G 1 )) 
a 

For a permutation x, define af to be the Card-Cyclic-to-Random shuffle starting at x. Our mixing 
time upper bound foliows from the foliowing lemma. 

Lemma 9 // permutations x and y differ by an adjacent transposition and n > A, there is a 
coupling of and a\ such that 

Ep(aZ,o-y)<e- a , 

where a = 2(log2 — log(e — 1)). 

Proof: There is another natural coupling of two Card-Cyclic-to-Random processes besides label 
coupling; we call this second coupling position coupling. In position coupling, the card is inserted 
into the same locations in both processes. Now assume that for some i < j, the permutation x can 
be obtained from y by transposing the cards with label i and j, as shown below. In the diagram, 
the kth X in the top row represents the same card as the fcth X in the bottom row. 

x: X X X i j X X X 
y: X X X j i X X X 

The coupling strategy is divided into 3 stages, corresponding to t in {1, ■ ■ ■ , i — 1}, {i, ■ ■ ■ ,j — 1}, 
and {j, ■ ■ ■ ,n} respectively. 

Stage 1: moving cards 1, . . . , i — 1. In this stage use position coupling. As is shown by diagram 
1 below, at the end of this stage we still have two permutations that differ only by a transposition 
of i and j. However, there may have been some cards inserted between cards i and j; we represent 
these cards with a's. 

af_! : X X i a a j X X 
crf_ 1 : X X j a a i X X 

diagram 1 

Stage 2: moving cards i, . . . ,j — 1. In this stage we use label coupling. At the end of this stage, 
some cards might have been inserted into the group of a's. We denote such cards with q"s. In 
addition, some cards might have been inserted between card j and the first X to the right of the 
card j. We represent them with 6's. Diagram 2 shows a typical pair of permutations after stage 2. 

a x j_ x : X X a a' a a' j b b X X 
ojf.! : X X j b b a a' a a' X X 

diagram 2 
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Stage 3: moving cards j, . . . ,n. Here we use label coupling again. Cards inserted into the group 
of a's and c/'s are represented with a**'s, and cards inserted into the group of 6's are represented 
with /3's. See diagram 3 below. 

<t£ : X X a a' a** a a** a' b ft b X X 
an ■ X X b f3' b a a' a** a a** a' X X 

diagram 3 

For t < n, let A t be the number of a's, a"s and a**'s, and let B t be the number of b's and /3"s, 
after card t has been moved. Note that 

Thus we are left to estimate E(A n B n ). 

Initially we have Aq = Bq = 0. Recall that in the first stage we use position coupling. For 
t < i — 1 we have B t = and A t satisfies 

F(A t+1 = A t \A t ) = n ' - 



n 



and 



This implies 



Hence 



F(A t+1 = A t + l\A t )-- ' " 1 



n 



E(A t+1 + l\A t ) = (A t + 1) ( 1 + 1 ) . (21) 



n 



n 



EAi-i = 1 + - - 1. (22) 



Recall that we use label coupling in the second stage. For i < t < j — 1, we have the following 
transition rule: 

F(A t+1 = A t ,B t+1 = Bt\A t ,Bt) = - - M - Bt - - , 

n 



and 



and 



This implies 



F(A t+1 =A t + l, B t+1 = Bt\A t , B t ) = 



F{A t+1 = A tl B t+1 = B t + l\A t , B t ) = h±l. 

n 



E(A t+l (Bt +l + l)\A t ,Bt) = At(B t + l) ( 1+ 2 



n 



Recall that B t = for all t < i — 1. Thus we have 



2 V~ l 

EA,-_i(Bj_i + 1) = EAi-i ( 1 + - J • (23) 
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Note that for t with i < t < j we have 



E(A t+ i\A t ) = At 1 + 



1 



n 



Thus EAj-i = EAj_i(l + i)-? -1 . Combining this with (23) and (22) gives 



EAj-iBj-i 



1 + 



n 



i-l 



1 + 



1 + 



For j <t < n we have the following transition probabilities: 

F(A t+1 = A t ,B t+1 = B t \A t ,B t ) 



n-At-Bt 



71 



F(A t+1 =At + l, B t+ i = Bt\A t ,B t 
F(A t+ i = A t , B t +i = B t + l\A t , B t 



n 
Bt 
n 



This implies 



E{A t+1 B t+1 \At,Bt) = A t Bt 1 + 



n 



Using (25), we obtain 



EA n B n 



1 + 



n 



i-l 



2 V"' 
1 + - 
n I 



1 + 



n 



1 + 



n-j+1 



(24) 



(25) 



Since 1 + ^ < (l + ^) 2 , the expression in square brackets is at most (l + % [\A + 3 ' — lj. 
Thus if we define /3 and 7 so that i = fin and j = jn, calculation yields that 

EA n B n < (eP - l)e 7 - /3 (e 7 " /3 - l)e 2(1 - 7) , 

if < p < log 2, and 

EA n B n < (e^ - l)e^(e^ - l)e 2 ^ (l + j ) , 

if log 2 < /3 < 1. The former expression is maximized, for 7 and /3 with < j3 < 7 < 1, by (^^) 2 - 
The maximum occurs when 7 = 1 and /3 = log^j-. Notice that log < log 2. Therefore, if 
a = 2(log2 - log(e - 1)), then 

E(A n B n ) < e~ a , 

for all < P < 7 < 1 and n > 4. which completes the proof. □ 



Proof of Theorem [2} We apply Theorem [8] to a round of the Card-Cyclic-to- Random shuffle. 
Since the diameter of S n with respect to adjacent transpositions is n ^" 2 - < n 2 , substituting the a 
of Lemma [9] into Theorem [8] gives 

*mix(e) < ; — 7, — r~7 7v(logn- 21oge). 

log 2 — log(e — 1) 



Acknowledgement. We thank Tonci Antunovic, Richard Pymar, Miklos Racz, Perla Sousi and 
Leonard Wong for many helpful comments. 



12 



References 

[1] R. Bubley and M. Dyer, Path Coupling: A technique for proving rapid mixing in Markov 
Chains, Proceedings of the 38th Annual Symposium on Foundation of Computer Science, pp. 
223-231, 1997. 

[2] P. Diaconis and M. Shahshahani, Generating a random permutation with random transposi- 
tions, Z. Wahr. verw. Gebiete, 57 159-179 (1981). 

[3] P. Diaconis and L. Saloff-Coste, Comparison techniques for random walks on finite groups, 
Ann. Probab. 21(1993) 2131-2156. 

[4] P. Diaconis and A. Ram, Analysis of Systematic Scan Metropolis Algorithm Using Iwahori- 
Hecke Algebra Techniques, Michigan Journal of Mathematics, 48(1):157-190. 

[5] W. Hoeffding, Probability inequalities for sums of bounded random variables. Journal of the 
American Statistical Association 58 (301): 13C30. 1963 

[6] D. Levin, Y. Peres and E. Wilmer, Markov Chains and mixing time, American Mathematical 
Society, Providence, RI, 2009. With a chapter by James G. Propp and David B. Wilson. 

[7] I. Mironov, (Not So) Random Shuffles of RC4, Proceedings of CRYPTO 2002, 304 C319. 

[8] E. Mossel, Y. Peres and A. Sinclair, Shuffling by semi-random transpositions, Proceedings of 
the 45th Annual IEEE Symposium on Foundations of Computer Science(FOCS'04) October 
17-19, 2004, Rome, Italy, 572-581, IEEE(2004). 

[9] R. Pinsky, Probabilistic and Combinatorial Aspects of the Card-Cyclic to Random Shuffle, 
preprint, 2011. 

[10] E. Subag, A Lower Bound for the Mixing Time of the Random-to-Random Insertions Shuffle, 
preprint, 2011. 

[11] L. Saloff-Coste and J. Zuniga, Convergence of some time inhomogeneous Markov chains via 
spectral techniques. Stochastic Processes and their Applications 117 (2007) 961-979. 

[12] J. Uyemura- Reyes, Random Walk, semi-direct products, and card shuffling, Ph.D. Thesis, 
Stanford University, 2002. 



13 



